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Abstract 

We have measured the intensity dependent rate and frequency shift of a photoassociation tran- 
sition in a quantum degenerate gas of "^Li. The rate increases hnearly with photoassociation laser 
intensity for low intensities, whereas saturation is observed at higher intensities. The measured 
rates and shifts agree reasonably well with theory within the estimated systematic uncertainties. 
Several theoretically predicted saturation mechanisms are discussed, but a theory in which satura- 
tion arises because of quantum mechanical unitarity agrees well with the data. 

PACS numbers: 03.75.Fi, 33.20.Kf, 33.70-w, 34.20. Cf 
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Photoassociation (PA) of ultracold atoms has been a remarkably useful tool for deter- 
mining scattering lengths characterizing ultracold atom collisions, for producing ultracold 
molecules, and for providing extremely precise measurements of atomic radiative lifetimes 
(see Refs. jl[ 0] for reviews). This utility is largely a consequence of the spectroscopic 
precision afforded by the small thermal broadening in a laser or evaporatively cooled gas. 
Quantum degenerate gases are especially interesting because the coherence of the atomic field 
may enable the formation of a molecular Bose-Einstein condensate (BEC) from an atomic 
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one by coherent Raman transitions [3|, y|, stimulated Raman adiabatic passage (STIRAP) 
0] , or by other coherent adiabatic population transfer schemes 0, 3 • 
There have been extensive theoretical studies on the rate of PA 

HIlllQ- The rate 

is predicted to increase linearly with intensity at low intensities, while various mechanisms 
have been proposed that cause saturation of the rate at higher intensities. Among these 
mechanisms are the quantum mechanical unitarity limit on the rate of atomic collisions 11[ , 
a break-down of the 2- mode approximation for PA of Bose condensates caused by cou- 



pling to non-condensed atomic modes 



and the depletion of the atomic pair 



correlation function 15]. Photoassociation resonances are also predicted to exhibit a spectral 
shift proportional to the light intensity caused by coupling to the continuum of free-atom 
states 17, 18, 3|- In contrast to theory, there are relatively few experimental measure- 
ments, and only two that could be considered "precise" (which we define as measurements 
with uncertainties of less than a factor of 2). We previously measured the spectral light 
shift using quantum degenerate ''Li and obtained good agreement with theory [2^. Both 
the spectral shift and the PA rate constant were recently measured in a Na condensate, and 
good agreement with two-body theory was found |21j]. Saturation was not observed in this 
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experiment. Saturation was observed in two other lower precision experiments 
these experiments were performed in a magneto-optical trap, where the temperatures were 
greater than 100 /iK and the corresponding unitarity-limited rates were quite small, of the 
order of 1 s~^ or less. 

We report precise measurements of both the rate of PA and the spectral shift in a quantum 
degenerate gas of ^Li atoms as a function of laser intensity. The rate is observed to saturate 
at the highest intensities, where the corresponding rate constants are nearly two orders of 
magnitude larger than any previously measured. The large PA rate is a consequence of 
quantum degeneracy, where the densities are relatively high and the unitarity limited rates 
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are large, as well as favorable free-bound state overlap in lithium. 

The apparatus will be discussed only briefly as it has been described in detail previously 
Atoms are confined by a magnetic trap with a bias field of ~1000 G. Atoms in the 
F = 2, mj7 = 2 hyperfine sub level of ^Li are cooled by evaporation to quantum degeneracy. 
Attractive interactions between atoms restrict the number of atoms in the Bose-Einstein 
condensate, Nq, to a relatively small number before the condensate collapses j^]. For the 
experimental conditions here, Nq ^ 1250, a small fraction of the total number of atoms. This 
fact, however, facilitates the achievement of repeatable temperatures, total atom numbers, 
and densities because quenching the gas below the BEG transition temperature Tc causes 
it to attain thermal equilibrium with a temperature T ~ Tc [26]. Therefore, following 
evaporation we allow the gas to equilibrate for approximately 12 s, at which point the total 
number of atoms ^ 8 x 10^, the gas is very near the BEG transition at T ~ 600 nK, and 
the peak density is ~4xl0^^ cm^. 

Following equilibration, the atoms are exposed to a pulse of photoassociation light. Up to 
85 mW of PA light, coupled out of a single- mode fiber is focused on to the atoms. The laser 
frequency is tuned to near resonance with the v' = 83 vibrational level of the l^Sj" state, 

n 

which has a binding energy of 60 GHz relative to the 2Pi/2 asymptote j27|]. The relative 
frequency of the PA laser is monitored by comparing it with a reference laser locked to an 
atomic transition using a scanning Fabry-Perot etalon. The v' = 83 level was chosen for its 
large ratio of the photoassociation rate to the rate of off-resonant scattering from the atomic 
resonance. 

Photoassociation can be detected because it causes a reduction in atom number when 
excited molecules spontaneously decay into pairs of hot atoms that escape the trap 
The number of atoms remaining following the PA pulse is determined by polarization phase- 



contrast imaging 



24| . and compared with a "background" measurement of the number 



of atoms without the PA pulse. Since imaging destroys the sample, the trap must be 
reloaded and the atoms evaporatively cooled for each image. Both and T are extracted 
by fitting the images to the distribution function for a Bose gas in an external potential. The 
photoassociation rate is determined for a given laser intensity / by adjusting the laser pulse 
length r to cause a 20% - 30% loss of the initial number of atoms. For the data presented 
here, r is between 12 /is at the highest intensity and 4 ms at the lowest. Other phenomena 
besides PA, such as spontaneous scattering from the atomic resonance and dipole forces, can 
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cause light- induced trap losses. Dipole forces, arising from gradients of the laser field, are 
minimized by using a large laser beam waist (1/e^ radius of between 220 and 320 fim). The 
background loss rate is measured by tuning the PA laser several GHz from the molecular 
resonance, and the data is adjusted accordingly, by up to 2% at the highest intensities. Fig. 
1 shows a typical photoassociation resonance curve for relatively low intensity. The curve 
is well- described by a Lorentzian with a width only slightly larger than the expected value 
of twice the atomic natural linewidth. Although we did not perform a systematic study of 
lineshapes, the resonances broaden as expected as / is increased. The resonance position 
shifts with laser intensity, as shown in Fig. 2. A linear fit to the data gives a slope of 
(-1.7 ±0.2) MHz/(W-cm-2). 

The measurement of the rate is expressed in terms of the on-resonance rate coefficient 
Kp, which is defined in terms of the time derivative of the density n{r,t) = —Kp{I)n'^{r,t). 
We neglect the effect of elastic collisions on the density distribution since r is much less 
than the characteristic time for elastic scattering of ~300 ms. The solution for the evolution 
of the density, n{r,t) = n(r, 0)/{l + Kp{I)n{r,0)t}, shows that the distribution becomes 
flatter during the PA pulse since the rate is largest at the highest densities The 
initial density n{r, 0) prior to the PA pulse is described by an equilibrium Bose-Einstein 
density distribution. The fractional trap loss, for a given value of Kp and r, is calculated 
by taking the ratio between the spatial integral of n{r, r) and the initial number of atoms. 
The best value for Kp is found numerically by matching the computed fractional loss with 
the measured one. Fig. 3 shows the dependence of the on-resonance rate coefficient Kp 
on the laser intensity /. In the low-intensity limit Kp is found to depend linearly on the 
laser intensity, with a slope of (7.9 ± 1.6) x 10^^*^ (cm'^-s""'^)/(W-cm^^). The rate constant 
saturates at high intensities, with a maximum value of (2.2 ± 0.2) x 10~^ cm^-s~^. The 
stated uncertainties are dominated by systematic effects. 

Several factors contribute to the uncertainties of the measured quantities. There is no 
direct way to know before exposure to the PA pulse and therefore the ability to measure 
PA-induced loss depends on the reproducibility of and T. We find that the statistical 
variations in A^ and T for the background measurements have standard deviations of only 
3% and 2%, respectively. Furthermore, the phase-space density, proportional to NT~^, is 
found to be within 5% of the critical value for all background measurements. Background 
measurements are typically made every fourth shot to monitor and make adjustments for 
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long-term drifts. Of more concern are systematic uncertainties. and T are subject to 
systematic uncertainties of 5% and 3%, respectively, due to uncertainties in the image mag- 
nification and probe laser polarization 2J]. Finally, the uncertainty in intensity is dominated 
by a systematic uncertainty in power measurement, which we estimate to be ~10%. These 
effects combine to give a 20% uncertainty to the slope of Kp, and a 10% uncertainty to both 
the rate constant maximum and the slope of the frequency shift. 



The data are compared with the theory of Ref. The transition dipole, dm-, for 

the free-bound transition to this high-lying vibrational level is simply related to the atomic 
dipole da- Since the binding energy of v' = 83 is large compared to the fine-structure splitting 
of the 2P atomic state and any hyperfine interactions, the electronic and nuclear spins of 
the atoms decouple from the PA light field. Therefore, for this H —>■ H transition, d^ lies 
along the internuclear axis. Averaging over the unknown orientation of the molecular axis 



relative to the polarization axis of the linearly polarized laser field gives a factor of \/l/3, so 



that dm = y'^/^da [28]. The free-bound Franck-Condon factors are numerically evaluated 
using molecular potentials determined previously Finally, because of quantum 

degeneracy, to calculate Kp we average the scattering matrix element over a Bose-Einstein 
energy distribution. 

The light-induced spectral shift is calculated to be —2.1 MHz/(W-cm~^), in reasonably 
good, but not perfect agreement with the measured value. The solid line in Fig. 3 shows 
the theoretical results for Kp. The dashed line is a low-intensity extrapolation of the theory 
which has a slope of 5.8 x 10^^*^ (cm^-s^^) / (W-cm~^), which again is in reasonable agreement. 
Finally, theory predicts saturation to Kp = 2.0 x 10^®cm^-s~^, which agrees very well with 
the measured value. 

The theory of Refs. [sl and show that Kp goes to a constant value, independent of 
T, in the low-T and low-/ limits. However, the free-bound matrix element depends on T in 
this theory. We find it convenient to define the free-bound matrix element as T-independent 
from the start, and use this to obtain a completely T-independent analytic approximation 
for Kp. Therefore, we scale the two-atom unbound wave function (pki'f') such that outside 
of the interior of the ground-state molecular potential </)fc(r) = ^^^^^Ip'^, where 6 is the 
s-wave scattering phase shift. The overlap between the unity-normalized molecular wave 
function and (f)k{f) is Mkv = 47r r'^(j)k{'r)ipl{r)dr. This matrix element has units 

of (length)^^'^ and reaches a constant value as T — ^ 0. In our experiment, for example, at 
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le classical outer turning 



29|, where Qq is the Bohr 



T ^ 600 nK and for v' = 83, the wave number ~ 2 x 10~^a~^, t 
point Rc = 103 tto and the s-wave scattering length a = —27.6 
radius. The value of |Mfc^,p for these conditions differs from the T ^ value by only 0.02%, 
showing that the system is well within the quantum threshold regime. At low intensities 
and temperatures, we find K„ = 77^|Mfc„P, where I sat and 7 are the molecular saturation 
intensity and spontaneous decay rate, respectively. For this transition I sat = 6/a, where 
la = 5.1 mW/cm^ is the atomic saturation intensity. 

Saturation comes about in the theory of Ref. llj because of the quantum mechanical 
unitarity limit on the rate of two-body collisions. Accordingly, the upper bound for the rate 
constant is = a(v), where a = -X^ is the maximum cross section between identical, 
non-condensed bosons, (v) is the thermally averaged velocity, and A is the thermal de Broglie 
wavelength. For T = 600 nK, iT^") = 2.5 x 10"^ cm^ s ^, in good agreement with the full 
heoretical calculation and experiment. The saturation observed in previous experiments 
2^ |2^ is consistent with the unitarity limit, although at rates many orders of magnitude 
lower than observed here. It is interesting to speculate on how our results relate to the 
other proposed saturation mechanisms. The reverse process of dissociation can result in 
the formation of pairs of "hot" atoms that are not returned to the original translational 
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state of the cold atomic gas [Ij, [IJ, [151, Hg. The energy width of the dissociated pairs is 
proportional to the PA rate, and results in a rate limit of r^hn^/^ /m Interestingly, at 
Tc, where A oc n '\ the unitarity limit gives the same maximum rate m. 
Since association must involve atom pairs in close proximity, depletion of the pair correlation 
function at short range may also limit the rate. From a classical, particle-like perspective, the 
maximum association rate is reached when atoms cannot move quickly enough to replenish 
the necessary close-range pairs, leading to a maximum rate of where d ~ n~^^^ is 

the mean separation between atoms and v is their mean velocity. But since v ~ h/mX and 
n ^ Nq/X^, the resulting limit is ~iVo ^^^hn"^^^ /m, which is again the same as the unitarity 
limit when T ~ Tc. The limit imposed by this classical arg ument, however, has already 
been violated in the BEC experiment reported in Ref. j2l|, in which a PA rate greatly 
in excess oi v/d was observed. Although no saturation was observed in that experiment, 
the rate was well below the unitarity limit, and Ur?!^ jm. Our experiment cannot clearly 
distinguish between the mechanisms because their predicted limits are all of the same order 
at Tc- However, the overall good quantitative agreement over a large range of intensity with 
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a theory whose rate is hmited by unitarity in the two-body scattering amphtude, strongly 
supports this explanation for our data. Finally, we point out that mechanisms limiting the 
PA rate will also limit the rate of molecule formation by magnetic Feshbach resonances. 

In summary, we report the results of precise measurements of the rate and spectral shift 
of photoassociation resonances of a quantum degenerate gas. The overall agreement between 
theory, with no free parameters, and our measurement is good. The agreement attests to 
the general validity of the theory and also points to quantum mechanical unitarity as the 
limitation of the rate for a gas at Tp. The distinction between the unitarity- limited model 
and the hot dissociation model could nevertheless be unambiguously tested with a lithium 
gas cooled significantly below Tc, an experiment we hope to perform in the future. 

We thank Robin Cote, Juha Javanainen and Paul Julienne for helpful discussions. This 
work was partially funded by grants from the National Science Foundation, the Office of 
Naval Research, the NASA, and the Welch Foundation. 
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FIG. 1: Photoassociation resonance curve for v' = 83. The solid circles are experimental data 
points and the solid line is a Lorentzian fit to them. The signal is the normalized number of atoms 
remaining in the trap after the photoassociation pulse of intensity 0.185 W/cm-^ and duration 2.0 
ms. The initial number of atoms is 5 x 10^ at a temperature of ~520 nK. The uncertainty in 
relative frequency is ~2 MHz, caused by drift of the PA and reference laser frequencies. The data 
fit to a rate constant with a Lorentzian width of 14 MHz, in good ag reement with the expected 
natural width of 11.7 MHz, which is twice the atomic natural width ll(J]. The thermal contribution 
to the lineshape is negligible at this temperature. 

FIG. 2: Light-induced frequency shift of the photoassociation resonance. The solid circles are data. 
The error bars are uncertainty estimates for the identification of the resonance peak positions. The 
dashed line is a fit to a straight line, where each data point is weighted by the inverse of the error 
bar length. 

FIG. 3: On-resonance photoassociation rate coefficient Kp as a function of laser intensity. The open 
circles are experimental data, obtained from trap loss analysis. The solid line is the theoretical 
prediction, with no adjustable parameters. The dashed line is an extrapolation of the theoretical 
result in the low- intensity limit. 
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Figure 1 of Prodan et al. 
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Figure 2 of Prodan et al. 
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Figure 3 of Prodan et al. 
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